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ABSTRACT 


Experimental data for both the charging and the blowdown process 
in a single air receiver is presented and interpreted in detail. The 
data indicate that heat transfer effects cause a radical departure 
from adiabatic behavior, and that such variances can be explained 
qualitatively on the basis of available simplified expressions for the 
State of a gas in a receiver. It is shown that the behavior of systems 
fulfilling the limiting conditions of these solutions can be adequately 
predicted and useful design results obtained. Methods for determining 
heat transfer convective conductances for use in the simplified solu- 
tions are discussed and evaluated, 

The experimental work was performed from January 1959 through 
April 1959 at the United States Naval Postgraduate School, Monterey, 


California 
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SYMBOLS AND UNITS 


English Letter §S ols 


Scale factor 

Area, £t2 

Width dimension, in. 

Constant 

Thermal capacitance of receiver shell, Btu/°R 
Specific heat at constant pressure, Btu/(1b°R) 
Specific heat at constant volume, Btu/(1b°R) 
Orifice diameter, in. 

Diameter of receiver, in. 

Mass velocity, lb/(hr £t2) 

Specific enthalpy, Btu/1lb 

Unit heat transfer convective conductance, Btu/(hr £t*°R) 
Jet momentum, defined in equation (18) 

Thermal conductivity, Btu/(hr ft? OR Et) 
Constant 

Length dimension, ft. 

Mass, lb. 

Pressure, lb/ft? 

Heat stored in insulation, Btu 

Heat transfer rate, Btu/hr 

Total heat transfer resistance, hr °R/Btu 
Universal gas constant/moleoular mass of gas, ft 1b/(1b°R) 
Time 

Absolute temperature, OR 


Total internal energy, Btu 





English Letter Symbols 


u Specific internal energy, Btu/1lb 
u Average velocity, ft/sec. 

V Volume, et? 

W Mass flow rate, lb/hr 

Ker Cylindrical coordinates 


Greek Letter Symbols 


oC Constant 
Temperature expansion factor, 1/°R 
Constant 
Constant 


Constant a 





Viscosity, lb/hr ft) 


SD =F BHF SY Hh 


Density, lb/ft? 

A Denotes a difference 
Nondimensional Grouping 

ce Pa C./ (Moc, ) 


a en) vahicn /c 


a 


M* « M/M, 


NTU #2 1/(R,; cy wo) = (hA),/c,, w 
NTU_« L/(R cy Wo) = (hA) /ey, Wo 


pe = P/P, 


= T/T, 


, dof 2 F 


= w/Wo 


Grashof Number 


q) 
| 
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( Lop * ghar/y°) 
viil 








Nondimensional Grouping 
Pr «= Prandtl Number (A c/ i) 


Re e« Reynolds Number COG) a) 


Subscripts 


o Refers to initial conditions 

1 Refers to inlet state 

c Refers to capacitance 

1 Refers to conditions inside receiver 

O Refers to environmental conditions outside receiver 
x Refers to conditions at a distance x 
0,0 Refers to conditions at the origin of a jet 
X,0O Refers to conditions at a point on a jet axis 
Xx, ¥ Refers to conditions at a point in an annulus 

m Refers to a mean value over a time pertod 
ave Refers to an average value at a given time 
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1. Introduction 

During blowdown and charging of a gas receiver, changes in the 
gas temperature may provide a substantial temperature difference for 
heat transfer between the gas and the receiver walls. In addition 
to heat transfer through the walls, the receiver mass itself may provide 
significant energy storage so that heat transfer from such a thermal 
mass becomes important. Since the thermal capacitance of the receiver 
is usually quite large relative to that of the contained gas, considerable 
heat transfer between the walls and the gas can occur with little 
change in the wall temperature. The usual approach in the design of 
charging and blowdown systems has been to assume an adiabatic process, 
Inasmuch as this procedure neglects heat traféfer effects, substantial 
error in the prediction of system behavior may result. 

When heat transfer rates are to be considered, it becomes necessary 
to specify the rate at which mass enters or leaves the receiver. Two 
types of flow of considerable interest are constant mass flow and flow 
through a critical flow nozzle. The former case finds considerable 
application in the blowdown wind tunnel. Here the problem consists of 
maintaining sufficiently constant stagnation conditions despite decreas- 
ing temperature from the expansion process in the reservoir. Methods 
applicable to this problem have been treated both analytically and 
experimentally by Murphy, et al., ou employing a separate thermal 
mass and air storage reservoir. Incorporated in these solutions are 
analytical expressions for specific modes of convection heat transfer 
between the air, reservoir walls, and thermal mass. Thus the value 


of these solutions is limited to systems behaving in a similar manner. 





lyumbers in brackets refer to references listed in the bibliography. 





Of perhaps more general interest are the solutions presented by 
Reynolds [2] for both constant and critical mass flows. These 
solutions employ a number of simplifying assumptions, the value of which 
are primarily dependent upon the accuracy of independently determining 
the heat transfer conductances between the gas, thermal mass, and 
environs. The assumption is made that the heat transfer resistances 
are uniform over all interior and exterior heat transfer surfaces and 
either invariant or independently predictable with time depending upon 
application. Critical flow blowdown experiments have been performed 
using time average free convection relationships which closely approxi- 
mate the conductances predicted from these solutions [3] . 

Little other published information is available for similar charg- 
ing and blowdown experiments. In receivers densely packed with heat 
capacitors, forced convection heat transfer correlations may apply. 

The fluid jet during charging may strongly influence the convection heat 
transfer mechanism. Imperfect mixing may occur during charging and 


cause temperature gradients within the system. 





2. Objectives 


Controlled single gas receiver charging and blowdown experiments 
were conducted at various mass flow rates and over varying ranges of 
thermal capacitance in order to accomplish the following objectives: 
(a) Evaluate previously developed simplified analytical expressions 
for the thermodynamic state of the gas as a function of time. 

(b) Determine the mechanism of heat transfer between the gas in the 
receiver and the receiver walls during these processes. 

(c) Investigate analytical methods of predicting heat transfer con- 
ductances for use with the expressions for the thermodynamic state of 
the gas. 

(d) Determine the feasibility of attaining an adiabatic process within 
a reasonable charging or blowdown time. 

(e) Investigate the temperature distribution in the receiver as 


affected by mixing during charging. 





3. Description of Test Apparatus 
3.1. General Description 


The test apparatus consisted of two 180 gallon pressure vessels of 
150 psig test, one of which contained an internal liner of 5/8 in. 
California redwood. The tanks were arranged as shown in Fig. 1 with 
charging air introduced at the bottom through various size orifice plates, 
Each tank was fitted at the top with alternative piping arrangements for 
either critical or constant mass flow blowdown as shown in Fig. 2, 
Physical dimensions are contained in Table 1, page 6 . Air for charging 
the tanks was supplied from an air bank through a pressure regulator 
set at 200 psig. Metal strips could be inserted vertically in the tanks 
to serve as thermal capacitors. Pressure and temperature were continu- 
ously recorded with time. A photographic view of the general system is 
contained in Fig. 3. 

3.2. Flow Metering 

Sharp edge orifice plates were available in increasing diameters 
ranging from 1/8 in, through 1/2 in. for interchangeable mounting in 
flanges attached to the tank top and bottom. Discharge coefficients for 
orifice sizes through 1/4 in, diameter were determined over a range of 
pressure ratios using a separate blowdown calibration system. These 
coefficients remained fairly constant up to the critical pressure ratio. 
Discharge coefficients for the larger orifices were not determined because 
of limitations in the calibration system. Instead, experimental data 
from actual charging and blowdown runs was used to determine average flow 
rates during constant flow. Constant flow charging runs up to 100 psig 
tank pressure were attainable from an air bank supply of 200 psig pressure, 


at close to ambient temperature, upstream from the orifice. Separate 
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arrangements permitted either critical or constant mass flow blowdown, 
The tank top could be fitted with a flange arrangement containing a 
seated orifice plate and a 1 in. quick opening gate valve for critical 
flow runs. Also available for constant mass flow blowdown was a l in. 
copper pipe and flange arrangement equipped with a hand valve and bourdon 
pressure gage upstream from the orifice plate for maintaining the desired 
pressure ratio. 

3.3. Temperature Measurements 

Each tank was equipped with six 30 gage copper-constantan thermo- 
couples. Four of these were arranged in series and spaced at equal 
volumes vertically in the tank to give an average tank air temperature, 
The two remaining junctions were available for monitoring wall or thermal 
capacitor temperature. Separate thermocouples were available for moni- 
toring both constant flow charging and constant flow blowdown air stag- 
nation temperature upstream from the orifice plates. Tank air and 
capacitor temperatures were recorded either on a continuously indicating 
Brown or on a Leeds Northrup potentiometer, 

3.4. Pressure Measurements 

Pressure was recorded from the output of a strain gage pressure 
pickup of the diaphragm type [};] . This pickup consisted of a seated 
diaphragm containing a center mounted SR4-A X 5 strain gage, together 
with two peripheral SR4-A5 gages as shown in Fig. 4. This arrangement 
was wired to form a four gage external bridge and provided a double out- 
put with temperature compensation. Calibrated sensitivity was 30 micro- 
inches per in. per psi. Output from the pickup was recorded on a Baldwin 
Strain recorder. A bourdon gage was mounted directly on the tank eee 
quick visual inspection of pressure. 


o) 





TABLE I 


Physical Dimensions of Experimental Apparatus 


1. Uninsulated Tank 

Volume: 23.9 ft? 
Internal area: 45.5 ft* 
Nominal internal diameter: 31 in. 
Nominal internal height: 54 in. 
Thermal mass: 835 lbs 

2. Insulated Tank 


Voliee. 17.0 £t- 


Internal area: 38,0 et? 
Nominal internal diameter: 29.5 in. 
Nominal internal height: 42 in, 

3. Redwood Insulation 
Thickness: 5/8 in. 
Density: 25 lbs/ft> 
Thermal Conductivity: 0.06 Btu/(hr ft °R) 
Specific heat: 0.6 Btu/(lb °R) 
Diffusivity: 0.004 f£t*/hr 

4. Strip Capacitors 
Material: 24 S Aluminum (Alclad) 
Mass: 10.0 lbs 
Area: 108 £t 


Gage: 0.012 in. 





4. Test Procedure 





Experimental work consisted of the following charging and blow- 
down runs: 

a. Charging and blowdown of the insulated tank under close to 
adiabatic conditions (no internal heat capacitors), but with some heat 
transfer to and from the redwood insulation, 

b. Charging and blowdown of the insulated tank with added heat 
capacitors and at constant mass flow rates to provide a finite capacitance 
with negligible inside heat transfer resistance. 

c. Charging and blowdown of the uninsulated tank at constant flow 
rates such that the tank walls provided an isothermal sink and source 
with a finite inside heat transfer resistance. 

In making a series of runs, the air vane was first charged with 
290 psig air from a reciprocating air compressor and allowed to cool to 
ambient temperature. The air bank was then drained of any condensed 
water, and the pressure regulator set at 200 psig. A quick opening valve 
between the tank orifice and pressure regulator was then opened for a 
constant mass flow charge through the selected orifice plate. The piping 
System between the air bank and tank was of sufficient thermal capacitance | 
that the air entered at nearly ambient temperature despite the blowdown 
process of the air bank. After reaching the desired tank pressure of 
about 100 psig, the tank was immediately blown down through either the 
critical flow quick opening valve or the constant flow hand valve arrange- 
ment. The latter process consisted of continually controlling the hand 
valve upstream of the blowdown orifice plate such that the upstream 
pressure remained constant at 30 psig. Once again, the piping and valve 
arrangement was of sufficient thermal capacitance that the temperature 
upstream from the orifice remained nearly ambient. 
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Throughout a typical cycle of charge and blowdown, tank air 
temperature was recorded at 5 sec, intervals with the Brown recorder 
and tank air pressure recorded continuously with the Baldwin recorder, 
Runs with added capacitors were repeated under identical conditions, 
recording the temperature of the capacitors with the Brown recorder, 
Runs varied in length from 11 to 240 seconds. Smooth curves were 
drawn through the recorded data plots. The flow rate was computed 
either from the usual critical flow metering equation using previously 
determined discharge coefficients or from the experimental mass-time 


data. 





5. Theoretical Considerations 





5.1. General 

The objective of analyses of blowdown and charging systems is to 
obtain relations expressing the thermodynamic state of the gas in the 
receiver as a function of time. It is usual to use temperature as the 
dependent variable with mass content as the independent variable. For 
a constant volume receiver, the pressure, temperature, and mass content 
are related by the perfect gas equation of state: 

PV «= M(R/m)T (1) 

Thus having found temperature as a function of mass the pressure may 
also be found as a function of mass, and knowing the mass-time relation- 
ship, all state functions may be expressed as functions of time. We 
shall be concerned in this investigation primarily with constant mass flow. 

General charging and blowdown solutions have been presented by 


Reynolds [2] using the following model: 


x ’ 
8 x 
x 
be > | amar alla 9 a 
x 
ts x 





These analyses assume that the state of the gas is uniform throughout 
the receiver and that the thermal capacitance of the receiver walls 

and other internal contents can be lumped into a single capacitance, 

Cc. = M. Cp The heat transfer resistance between the capacitance 

and the environment is denoted by R,, and the resistance between the 
capacitance and the air in the tank by Ry. The mass of gas in the tank 


is M and the mass flow rate is w. 
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Assuming that R,, , Ri» and C. are invariant in time, a combination 
of energy balance and heat transfer rate equations leads to the general 


equation! for charging and blowdown at constant mass flow of the form: 
2 
ad d BG apes * * 
M > +¢24, +7" dT * , 
“%2 * MM a + TT + e Pa O ( ) 
dM dM d M* 0 
The constants | ro Rian and @ are functions of the system parameters 


NTU, NTU.) , and C*.. These parameters are defined as follows: 


i eee ee CRA (3) 
i ov Wo Wo Cy 
ce a) Cy ms Uy Cy 

on = Me ¢p = Ge (5) 
M, Cy Mo Cy 


The NTU groupings are the familiar "number of transfer units" used 
in heat exchanger analysis and represent dimensionless conductances 
(reciprocal of resistance) and also serve as a measure of the rate 
of heat transfer of the process. The Ce, parameter represents the ratio 
of the capacitance of the receiver walls and internal thermal mass to 
the initial capacitance of the gas. 

While equation (2) may be solved by analytical, graphical, and analog 
pomputer- methods, the use of such solutions for engineering applica- 
tions is often difficult or impracticable. In the majority of engineering 
applications the magnitude of the system parameters NTU, NTU O°? and 


C,* may be such that more useful closed form solutions may be obtained. 





1 See Appendix I for the derivation of the general equation. 


2 An analog computer method for the solution of the general equation 
for blowdown is presented and evaluated in Appendix III. 
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Several of these solutions have been developed by Reynolds [2] 
from basic principles, and it is with the application and evalua- 
tion of these solutions for constant flow that this thesis is 
concerned. Only the resultant expressions for temperature as a 
function of mass content will be given in this section for the 
application of interest. For convenience, the derivations of the 
general and of the simplified solutions are contained in Appendix I. 


2.2. Adiabatic Charging and Blowdown 


If there is no heat transfer to or from the gas, the process 


wr 


within the tank is adiabatic. An adiabatic process might be approxi- 
mated during a very rapid charge or discharge, or during a relatively 
Slow process with a receiver lined with an insulator between the gas 
and wall capacitance. Adiabatic behavior corresponds to the situation 
where C. * = O and NTU ,, = 0, or when NTU = 0. Since no heat transfer 
is involved, the system is independent of mass flow rate. The follow- 


ing simple model and solutions apply: 


ive <> w 
charging % 
qe = hap = Anat (6) 


b lowdown ae i ek D (7) 
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5.3. Charging and Blowdown with Negligible Inside Resistance 


Systems having a high inside conductance (negligible inside heat 
transfer resistance) combined with a finite but low capacitance are 
often found in service. The limiting condition of NTU = cO - implies 
that the capacitance temperature is identical with the gas tempera- 
ture. This might be the case with the charging or blowdown of a thin 
walled pressure vessel where jet flow and gas currents create con- 
siderable turbulence resulting in high inside conductance values. The 


following simple model and solutions apply: 





(<3 fn Sen (8) 
we Byes C+NTUo- kT — NTUOTRIACZ +MY tT" +NTUy To 
I+ NTU 
b lowdown Ce ye a 
re (K-14 NTUo-Te NTU) \_ 14+ Co / + To NTUg 
K-|+NTU, 


One real value of these solutions is that they supply definite infor- 
mation for determining the effects of capacitance on the behavior of 
a system. 

9.4. Charging and Blowdown with an Isothermal Sink and Source 

In many systems the thermal capacitance of the metal receiver 
walls far exceeds the thermal capacitance of the gas. In such a case 


the temperature change of the walls is much less than that of the gas, 
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and in the limiting case of c.™ = o , heat transfer occurs only 
between an isothermal sink or source and the gas. This might ordinarily 
be the case with any relatively thick walled receiver or a receiver 
packed with thermal capacitors. The following simple model and 


solutions apply: 





Wes 
Wo 
charging nA ia 
=", KTPANTU qo — OKT le ntu entu Tc © 
_ }+NTU 
blowdown kK-I4+NTU * 
*. K-!+NTU-Te NTU) M™ fuilc NIU Ge 
- kK-1+NTU : 


The value of these solutions is further enhanced in that they can be 
arranged to give stepwise changes in T* when the parameters NTU and 


Te are varying with time. 


charging Pa a NTU 
AT = KT 4NTU Te ales |-(Mr OM) (4 (12) 
J+ NTU mM* 
blowdown x\ K-It+ NTU 
x | NTU M’+AmM _ 13 
atte tts Te (eau || 3 


5.5. Isothermal Case 


Receivers having a large thermal capacitance accompanied by high 
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values of inside heat transfer conductance may exhibit essentially 
isothermal behavior. This might be the case with a blowdown wind 
tunnel filled with a metal matrix for the purpose of obtaining a 
near-isothermal blowdown. The effect is to make both C,* and NTU 


very large. The following simple model applies: 


fa, 

alias. TM << on  e 
é 

a4 4 


44 4 
raed 


The temperature ratio for both charging and blowdown is by definition: 


T* = 1 (14) 


5.6. Application 

While it is not true that all engineering systems involving charging 
or blowdown processes may be treated with these special solutions, 
many of eee closely fit the limiting conditions. To obtain useful ’ 
results with these constant resistance solutions, one must use a suitable 
average NIU. In any real receiver the wide variations in flow rates, 
temperatures, and pressures will certainly cause considerable variation 
in heat transfer convective conductance. The remainder of the report 
will treat in detail the problem of obtaining suitable values of average 
or time varying conductances for applying these solutions to typical 


systems. Theoretical considerations of heat transfer conductances will 


be discussed where directly applicable to the results. 
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6. Summary of Experimental Results 


6.1. Form of Results 

Data for all the charging and blowdown runs consisted of tank 
air temperature and pressure measurements as a function of time. From 
these data, P¥ and T* were determined from definition and M* computed 
from the relationship, P¥ = M* T* obtained from equation (1). All 
graphical results are in the form of T* vs. M* since T* is more sensi- 
tive than P* to the effects of heat transfer. Smooth curves were 
drawn through the plotted data with the comparable theoretical solution 
indicated by a dashed line, The adiabatic solutions are also shown 
since any departures from these solutions are an indication of the 
occurrence of heat transfer. 


6.2. Uninsulated Tank 





A series of constant flow charging and blowdown runs were con- 
ducted with the uninsulated tank in the "as found" condition with no 
additional heat capacitors. The capacitance of the tank walls was 
of sufficient magnitude that heat transfer essentially occurred to 


and from an isothermal sink and source. The following runs were 


conducted; 
TABLE 2 
Summary of Uninsulated Tank Runs 

Run Orifice Wo t cx 

No. Type Dia. (in) (lbas/hr) (sec) 9 NIU 
l Charge 1/8 180 200 300 10.5 
2 Charge 3/16 380 100 300 7.95 
3 Charge 1/4 711 55 300 5.7 
4 Charge 5/16 1140 35 300 al 
5 Blow 3/16 105 240 41 5.0 
6 Blow 1/4 166 140 41 ie 
7 ~~ Blow 5/16 257 80 41 2.9 
8 Blow 3/8 375 70 41 2.6 
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The test results are presented in Tables 11 and 12 of Appendix IV 
and are shown graphically on Figs. 5, 6, and 7. Solutions based upon 


equations (10) and (11) to fit the experimental results are shown graph- 
ically for comparison. The values of NTU listed in Table 2, page 15, 
are those used in these equations. The values of Ce were obtained direct- 
ly from equation (5) based on a wall thermal capacitance of C, = 92 Btu/R. 
6.3. Insulated Tank 

Several charging and blowdown runs were made with the insulated 
tank without added heat capacitors in order to attain, as nearly as 
possible, an adiabatic process and to gage the effectiveness of the 


insulation. The following runs were conducted: 


TABLE 3 
Summary of Insulated Tank Runs 
Run Orifice Wo t 
No. Type Dia. (in (1b/hr ) gec 
9 Charge 5/16 1130 20 
10 Charge 3/8 1320 16 
ll Charge le 1960 ul 
2 Blow 1/4 395 120 
13 Blow 3/8 844 LO 
14 Blow 7/16 860 35 
iS Blow We 1100 26 


The test results are presented in Tables 13 and 14 of Appendix IV and 
are shown graphically on Figs. 8 and 9. The adiabatic solutions from 
equations (6) and (7) are shown for comparison. 
6.4. Insulated Tank with Added Heat Capacitors 
A series of constant mass flow charging and blowdown runs were made 


with a set of vertical strip aluminum heat capaaitors inserted in the 
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insulated tank. Capacitor configuration and air flow rates were 
selected such that the air temperature closely followed that of the 
capacitance for the case of heat transfer with negligible inside 
resistance. The following runs were conducted: 


TABLE 4 


Summary of Insulated Tank Runs 
with Added Capacitors 





Run Orifice Wo 2 

No. Type Dia. (in (ib/hr) = (sec) O 
16 Charge 3/16 4O7 70 10.0 
17 Charge 1/4 704 40 10.0 
18 Charge 5/16 1104 25 10.0 
19 Blow 5/16 255 70 1.44 
20 Blow 3/8 356 50 1.44 
21 Blow wile 497 35 1.44 


The test results are presented in Tables 15 and 16 of Appendix IV and 
are shown graphically on Figs. 10 and 11. Solutions based on equations 
(8) and (9) to fit the experimental results are shown ePephilcally for 
comparison. The values of or were obtained direotly from equation (5) 


based on a thermal capacitance of Co = 2.2 Btu/OR. 
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7. Discussion of Results 

7.1. Uninsulated Tank 

7.1.1. General 

Several charging and blowdown runs, as listed in Table2, page 15, 
were made with the uninsulated tank in the "as found" condition with 
no additional heat transfer surfaces. The tank was selected as a 
typical medium pressure air receiver of the type utilized for com- 
pressed air systems, The cylindrical tank walls and dished heads 
.were constructed of 1/2-in. steel plate having a total thermal capa- 
citance of C, = 92.0 BIU/°R. 

7.1.2. @harging of Uninsulated Tank 

The mass of air in the tank prior to charging was M, = 1.8 lbs. 
for all runs, giving a value of C.* = 300 from equation (5). This 
value is well above the figure of c.* « 40 recommended by Reynolds [2] 
as the minimum for infinite capacitance behavior. Therefore, the 
solution for charging with heat transfer to an isothermal sink would 
be expected to apply. The capacitance temperature during all charging 
runs remained essentially constant at 533 °R giving a value of 
(hs = 0.995. Charging air entered at a constant mass flow rate at 
a temperature of 531 °R giving a value of T,* = 0.990. Inserting 
these constant values in equation (10) together with the best value 
of NTU fitting the experimental points, resulted in the solutions of 
T* vs. M* shown graphically on Fig. 5. For values of Ma> 3 it was 
found that equation (10) yields a constant value of T*. Inspection 
of the experimental data shows that for all runs a constant value of 


T* was attained after about 35 seconds. The valves of NTU used in 
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equation (10) were therefore selected to give a calculated T* coincid= 
ing with these constant experimental values. 


A convective heat transfer conductance, h, can be found from the 


(hA)i 


UI, Cy 2 


definition: NTU = 


If the best single value of NTU fitting the experimental points is sub- 
stituted in this definition, the resulting value of conductance, h, 
becomes a mean convective conductance, h,, implied constant over the 
entire charging run. Values of h, from NTU are listed in Table 5 below. 


TABLE 


Qharging Qonductances 
Uninsulated Tank 


MAX he MAX. hove 
SE ales atl ASSEI (erncitn, euinly 
1 180 10.5 Te2 261 fel 
2 Ze 105 10.7 2.5 11.0 
5 711 ore 15.0 2e/ 14.0 
4 1140 4A aes: 209 16.5 


If the mechanism of heat transfer between the air in the tank and the 


tank walls is simple free convection, McAdams [5] recommends as a 
correlation 


Neo = “ 0.13 (Gr P,)” (15) 


10 11 
for 10 “<Gr Pr <10 . Inserting in equation (15) the maximum value 
of temperature attained during each run results in the maximum values 


of average convective heat transfer conductance, h listed in Table 5. 


ave’ 


Since these conductances are considerably lower than those calculated 


from NTU, it is therefore concluded that the heat transfer mechanism for 
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charging is forced convection. It is also apparent that the conduc- 
tances vary with time since the experimental values of T* decrease 
after attaining an initial peak value. This bears out the observation 
of Reynolds [2] that in any real receiver the wide variations in 
temperature and pressure would cause considerable variation in the heat 
transfer convective conductance. 

The presence of forced convection considerably conpounaemtHe 
problem of predicting heat transfer conductances. There is inevitably 
present a jet of fluid discharging through the orifice into a region 
of increasing pressure. While no analytical treatment of this specific 
problem has been found, there is considerable information available 
concerning the behavior of free jets discharging into a region of 
fluid at rest, Theory concerning the transport of momentum in an iso- 
thermal, turbulent, free jet of air discharging at subsonic velocities 
into air at the same density might be expected to lead to useful results 
in the present problem, 

Therefore, let us consider the following to be a simplified model 


of the jet and resulting flow inside the tank: 





Considering the jet as originating from a point at the orifice, obser- 


vations show that such a jet spreads out conically with a width b 
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directly proportional to the distance x from the origin [6]. The 
cone angle has been found to be about 25° [7] independent of the veloc- 


ity at the origin, u As the jet advances, the maximum velocity, 


ono ii 
Ux 0? at a point x on the centerline diminishes proportional to 1/x [a]. 
Fresh masses of fluid are continually being drawn in so that the mass 
flow at succeeding cross sections is not the same. The jet momentum, 
however, is constant since the pressure is assumed the same as that of 


the surrounding fluid fs]. These relationships may be expressed as 


follows: 


bm Xx (16) 

tro (17) 
ae, HE: 2 

| cok ae OSL Oe 


These relationships have been combined and experimentally verified 
by Alexander, et al., [9] in the form of momentum flux ratios expressible 
in terms of the ratio of orifice diameter to the distance x, as follows: 

(Zia =(a)" (19) 
(Pa ews K x 
where K is a constant experimentally determined for various initial 
velocities. Since the orifice discharge region pressure is always less 
than the critical pressure, the air will have expanded at the orifice 
throat only to the critical pressure. Therefore the quantity (PT)o, 6 
will be a constant independent of orifice size during constant flow 


charging. Equation (19) may be placed in the form, 
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= HO) oe J oe C 





(20) 


Where is now dependent only on the tank pressure and temperature, 
assumed oonstant throughout the tank. From the familiar cfitical flow 
metering equation and from continuity, an average discharge velocity 

was found to be 1050 ft/sec. At this velocity, Alexander [9] recommends 


K20.134, giving the following values of C from equation (20), 


a Nai c (ft 1b) 1/2 /sec 


1/8 750 
3/16 1125 
1/4 1500 
5/16 1690 


from which Lx 0 may be found at any point along the jet axis. 

Now consider the jet as impinging upon the dished head of the tank 
and being deflected, as shown in the model, down along the walls in a 
concentric annulus. If this be the case, then in order that there be no 
pressure difference throughout the tank, the mass flow upwards in the jet 
must equal the mass flow down through the outer annulus at any cross 
section. From continuity, the average velocity in the outer annulus at 


any distance x in terms of the annular and jet areas is 


U = ie UW 
,O 
x,/L DS ore x (21) 


For turbulent flow in concentric annuli, McAdams [ 5 | recommends for the 


local convective heat transfer conductance, 


hy = 023 (22) 


fa ak - 
(Pr)? (Re) 


where Re is based on the equivalent diameter, 
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D, = D~b, (23) 


Substituting equations (20) and (21) in equation (22) and using average 


values for the Prandtl number gives 


4 ie 
hy = 26 ge, Gane, el 


- VD, A «A ~ em 


Equation (22a) was used to calculate the local convective cond- 


Se. 
etantagn re D ) Sret en cP  h t  h ORETA 


Nene 





uctance, hy» at equally spaced intervals along the walls at various 


sn nomenon he 
et cneatnee er 


time intervals for each charging run using the experimental pressure- 
temperature data to evaluate the density. Thus for each run a series 


of curves of the type shown below were plotted. These curves show a 


h, (tu/nft? °R) 


O Se er L 


decreasing local conductance, tank top to bottom. From these curves, an 
average convective heat transfer conductance, Neve? was determined for 
each time increment. Values of have are shown graphically on Fig.(12). 

In all cases the predicted hes was found to have occurred at a point 
about $ the total tank height. For comparison, equation (10) was solved 
for an NTU corresponding to each experimental T ve. M* point. Using these 
values of NTU, an experimental Ney was calculated from the definition 


@ 


of NTU for each experimental point. These experimental values of Ms 


are also shown on Fig. (12). 
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Comparison of the curves on Fig. (12) shows a good correlation 


between the predicted h,,, and that determined from the experimental 


e 
points. As expected, increases in pressure and density predominate at the 
higher flow rates giving increased convective conductances and an increased 


NTU. The predicted values of h,,, follow the trend of the experimental 


é 
values extremely well although lagging somewhat in magnitude. At the 

lower flow rates the average convective conductance remains almost constant 
and this trend is not predicted as well by equation (22a). At such flow 
rates, however, the process is essentially isothermal and therefore the 
simple relationship T* = 1 gives adequate results. At the higher flow 
rates where NTU is changing with time equation (22a) can be used to predict 
NTU for use in stepwise equation (12) with results at least adequate for 


design and certainly better than would be predicted from adiabatic behavior. 


7.1.3- Blowdown of Uninsulated Tank 





The mass of air in the tank prior to blewdown was M, « 13.2 lbs. 
for all runs yielding a value of C,* = 41 from equation (5). This is well 
above the value of C, * = 6.5 recommended by Reynolds [2] as a minimum 
for infinite capacitance blowdown, Therefore, the solution for blowdown 
with heat transfer from an isothermal source would be expected to apply. 
The capacitance temperature during all blowdown runs remained essentially 
constant at 535 °R, Solutions of T* vs. M* based on equation (11) are 
shown graphically on Fig. 6 and 7 for the best value of NTU fitting the 
experimental points. 

Once again, solving for the mean heat transfer conductance over the 
entire run, h,, from the definition of NTU yields the values listed in 


Table 6, 
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TABLE 6 


Blowdown Oonductances 7 
Uninsulated Tank 


he 


(from Pe (free ott) 


Run W NTU % 
No. (lb/hr) (Eq. 11) (Btu/hrft“°R) (Btu/hrft“°R) Difference 


5 105 5.0 1.96 1.80 8 
6 166 3.5 2016 1.85 14 
7 257 2.9 2.79 1.90 32 
8 i 2.6 3.65 1.90 hg 


It is now of interest to compare the test results with conductances 
predicted on the basis of an assumption of a steady state turbulent free 
convection boundary layer over the entire heating surface. If average 
properties of air are inserted in equation (15), but the pressure depend- 
ence of the density retained, we obtain the equation for an average 


convective conductance, 
Vi 
hiade = 028 Ceaiae ss (15a) 


where P is the pressure in psia and AT is the temperature difference 
between the air and the tank walls in °F. Representative results from 
this equation for runs 5 and 6 are plotted on Fig. 13. In all cases the 
conductance drops sharply to zero when the air temperature reaches that 
of the tank wall and then rises rapidly to remain relatively constant 
over the remainder of the run. The mean predicted conductances with time 
from these curves are shown in Table 6 above. 
In each case, the predicted mean conductances are less than the 

measured mean conductances. At the slower flow rates the per cent dif- 


ference is reasonably small considering the inherent uncertainty in 
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matching the analytical curves and the experimental data. Reynolds [3] 
found differences of the same order of magnitude for similar blowdowns 
at critical flow, Runs 7 and 8, however, show that this difference 
becomes more pronounced as the flow rate increases. Steady state 
conditions still seem to persist since the experimental points for 
these runs (Fig. 7) follow the theoretical curves down to a point where 
constant flow departure begins. Large flow rates, however, yield small 
values of NIU and a closer approach to adiabatic behavior (Fig. 7) so 
that the precise evaluation of conductance might not be required. If, 
as Reynolds Ey points out, it is the relatively slow blowdown process 
that provides the most diffiaulty for the engineer because of heat transfer, 
then for such cases steady-state heat transfer correlations apparently 
provide a good approximation. | 

Blowdown systems where the evaluation of forced convection con- 
ductances becomes important are common enough, Murphy (1] found that 
for a blowdown wind tunnel with added heat capacitors, heat transfer 
occurred at a rate of about eight times that predicted for free convec- 
tion. This system differed from the simple case considered here in 
that capacitors were installed in a section separate from the air receiver, 
through which the air passed during blowdown and charging. Evaluation 
of forced convection conductance, therefore, would be dependent upon 
design characteristics and presumably predictable from available theory. 

For a tank of dimensions similar to the one used in these experi- 
ments, it is probable that the blowdown heat transfer mechanism would 
be free convection. Thus the relationships similar to equation (15) 


would give good approximations especially if increased slightly by 
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good judgment to allow for some forced convection effects at higher 
flow rates. 
7.2. Insulated Tank 

7.2.1. General 

Several charging and blowdown runs, as listed in Table 3, page 16, 
were conducted with the insulated tank in order to determine how closely 
the adiabatic solutions would apply to a system deliberately insulated 
from the containing walls. Choice of an insulation material presented 
a compromise of several considerations. Ideally, a material with both 
a low diffusivity (k/ cy) and low thermal storage capacity (Pe,) 
is required in order to have the surface temperature follow that of 
the air while the change of energy storage remains small. Inasmuch as 
the tank was relatively large in size, a material of sufficient struc- 
tural rigidity was also needed for ease of installation and dimensional 
stability under pressure. Foamed polystyrene insulation (Dow Chemical 
Styrofoam, HD2) was first selected and installed in the tank. This 
material, although possessing both a low diffusivity and low energy 
storage capacity, buckled under a hydrostatic pressure of 50 psig and 
came loose from the tank walls. This relatively low pressure also 
reduced the thickness of the Styrofoam about 10 percent. Therefore, 
redwood was selected as the most economical alternative material despite 
a relatively high density (25 lbs/ft’) and resulting high thermal 
storage capacity. 

7.2.2. Charging of Insulated Tank 

Reynolds [2] predicts that the adiabatic charging solution should 


be applicable for values of NTU < 0.25 for all magnitudes of C\*. From 
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equation (3) it is seen that a low NTU for a given area is attained 
by the combination of a high initial flow rate and a low heat transfer 
conductance. Therefore, orifice plates were selected to give the 
highest flow rate possible within a reasonable charging time. The 
results of the three runs, as displayed on Fig. 8,show considerable 
departure from the adiabatic solution. 

Inasmuch as the flow rates and areas are fixed, let us consider 
the magnitude of the mean convective heat transfer conductance, h, 
required to give a maximum NTU = 0.25. Solving for h, from the 


definition of NIU yields the conductance values listed in Table 7. 


TABLE 7 


Charging Conductances - Insulated Tank 





Run + MAX. "p 
(for NTU = 0.25) 
No. (1b7hr) (Btu/hr ft 20R) 
9 1130 1.26 
10 1320 1.48 
11 1960 "2,19 


It was seen, however, from the uninsulated tank results of Table 5, 
pe. 19 , that the conductances are in reality an order of magnitude 
greater than the above maximum values. This would then give a 
corresponding value of NTU w= 2.5. Even though Fig. 8 shows a nearer 
approach to the adiabatic line for the fastest run, it is conceivable 
that at faster charges the heat transfer conductances would increase 
to such an extent as to cause a departure from rather than a con- 
tinual approach to the line. 

Reference [2] also gives ag a criterion for adiabatic charging 
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the condition of ©. <0.04 for NTU reasonably low. Even if 
NTU # 2.5 is reasonably small, it will be shown in the next section 
that the G. of the insulation increases with time to values consider- 
ably in excess of the maximum Q” = 0.04, This low value of Co is 
unrealistic at best since its attainment would require a capacitance 
in this case of only C, = 0.086 Btu/OR ( equation (5) with Mo » 12.8 lb), 
a low value for the best of insulation. 

Thus it would seem that the criteria of NTU <0125, al though 
well within the range of free convection heat transfer conductances, 
is unattainable with a tank where the charging jet effect raises the 
conductance materially. As a result, it was not possible to attain 
a small ron are while charging the insulated tank, to evaluate 
fully the limiting conditions for adiabatic charging. 

aces Blowdown of Insulated Tank 

For the case of adiabatic blowdown, Reynolds [2] limits the con- 
ditions to NTU <.0.08 for all magnitudes of oe As in the case of 
charging, small values of NTU were sought by increasing the flow rate. 
Therefore, all blowdown runs were made at critical mass flow in order 
to attain the greatest value of initial flow rate. These runs are sun- 
marized in Table 3, page 16, and shown graphically on Fig. 9. As in 
the case of charging, Fig. 9 shows a considerable departure from adia- 
batic behavior. If once again we solve for the mean convective heat 
transfer conductance, h,, from the limiting value of NTU = 0.08, we 
obtain the values listed in Table 8, page 30. 

Results previously obtained from blowdown runs conducted with the 


insulated tank (Table 6, page 25) indicate that conductances slightly 
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TABLE 8 


Blowdown Conductances 
Insulated Tank 


Run Wo MAX Ams 
Now (1bs/hr ) (Btu/hr ft °°R) 

12 395 0.14 

Ue 841 0.30 

14 860 0.31 

15) 1100 0.29 


greater than those predicted for free convection might be expected. 
Inasmuch as these free convection conductances are an order of magnie- 
tude greater than the values of h, listed in Table &, the actual NTU 
for blowdown of the insulated tank is considerably greater than the 
value of NTU = 0.08. 

Several phenomena of interest are deducible from Fig. 9. It is 
noted that run 12, although at a slower flow rate and presumably at a 
higher NTU than the other runs, approaches the adiabatic line early in 
the run. This Micht be explained qualitatively by considering the test 
procedure. Blowdown was started immediately after charging and thus the 
air was at a relatively high temperature while the insulation remained 
essentially at ambient. Therefore, for a long blowdown there was time 
for considerable heat transfer to such an isothermal sink, dropping the 
temperature appreciably. At the end of the run conditions were reversed. 
The insulation now acted as an isothermal source for heat transfer from 
the lower temperature air. Thus toward the end of the run the air 
temperature dropped less rapidly. (This effect is also quite apparent 


with blowdown runs 5 and 6, Fig.6, conducted with the uninsulated tank.) 
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As flow rates increase and the run time decreases, this ‘apparent 
adiabatic" effect becomes less pronounced and, in the case of run 13, 
the temperature drop appears more nearly linear and departs further 
from adiabatic. As the flow rate increases further this linearity is 
more apparent and, since NIU decreases, the behavior again approaches 
adiabatic as seen by run 14. A still further increase in flow rate 
as with run 15 might be expected to further approach the adiabatic line, 
but this was not the case. Here, it is probable that h was also increas- 
ing due to increased convection, thus offsetting the increase in w, and 
in effect increasing NIU. Thus it would appear that a limiting approach 
to adiabatic behavior is all that may be attained within reasonable 
blowdown times with an arrangement such as this. 

Again, R@ynolds [2] establishes the condition of c.* ¢ 0.04 and 
NTU reasonably small as a criterion for adiabatic behavior. These con- 
ditions might conceivably become more realistic for blowdown than for 
charging since the initial mass of air in the tank is considerably 
greater than for charging and NIU small in comparison. In this case, 
an M = 8.8 lbs. yields a C, = 0.06 Btu/°R, a small value considering 
the capacity of the redwood insulation for storage of energy. 

7.3. Added Heat Capacitors in Insulated Tank 

7.3.1. General 

Several charging and blowdown runs, as listed in Table 4, page 1/7, 
were conducted with added heat capacitors inserted in the insulated 
tank. The insulated tank was used in order to eliminate heat transfer 
from outside the receiver (NTU LY = 0) and to minimize the effects of 


the tank walls so that most of the heating surface would be the heat 
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capacitors. Sixteén 0.012 in. thick aluminum strip capacitors, each 
12 ine by 40 in.w, were installed in a vertical and radial position in 
the tank and spaced evenly in order to interfere as little as possible 
with the air flow patterns previously discussed. The surface area of 
the capacitors was 108 ft? inoreasing the total area, including the tank 
walls, to 146 rt®, Total weight of the capacitors was 10.0 lbs. Thus 
the primary effect of the capacitance was to increase the heat transfer 
area and, therefore, NTU. 
7e3e20 Charging of Insulated Tank with Added Capacitors 

The mass of air in the tank prior to charging was M, = 1.3 lbae, 
yielding a value of C* «10.0 from equation (5)e Inasmuch as this value 
ou o” ia considerably lower than the minimum of oR = 40 recommended for 
infinite capacitance behavior [ele it is necessary now to estimate NTU. 
Since an approximation is justifiable, the conductance results prev- 
iously obtained with the uninsulated tank were used. Insertion of the 
minimum experimental values of average conductance from Fig. 12 in the 


definition of NTU yields the values of NTU listed in Table 9 below. 


TABLE 9 


Charging Conductances 
Added Capacitors in Insulated Tank 


mine have . 
Run Wo (Fig. 12) NTU Go 
No. (lb/hr) (Btu/hrft2°R) (Eq. 3) 
16 407, 8.0 16.9 10.0 
17 714 10.3 12.6 10.0 
18 1104 12.0 9.4 10,0 


These values are each above the minimum value of NTU> 7 recommended by 


Reynolds [2] for infinite NTU behavior. Inasmuch as the temperature of 
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the capacitance closely followed that of the air, the selution for 
negligible inside resistance would be expected to apply. 
If NTU «= O and Ti = 1, equation (8) reduces to: 


a 


ap 14-04 (ee ar] (8a ) 

Equation (8a) solved for O° #10.0 is shown graphically on Fig. 10. 
Comparison of this curve with the experimental mean plot shows a cone 
siderable variance. Inasmuch es equation (8a) is independent of a heat 
transfer rate (i.e. NTU,* O), it is also independent of mass flow rate 
and all runs should follow the same curve. There is some variance here 
as seen in Fig. 10. 

Since equation (8a) is sensitive to small changes in Co, any in- 
crease in receiver capacitance, C,, with time would considerably affect 
tT. This is especially true since the initial capacitance of the air 
(Mo Cy ) is relatively small. As was seen from the insulated tank runs 
without heat capacitors (Fig. 5), there is considerable heat transfer to 
the insulation during a typical charge. Therefore, a comparison of the 
increase in insulation of! alone during a run with that given by equation 
(8a) for the experimental data should give an indication of the applio- 
ability of this solution. 

The temperature distribution in the redwood as a function of time 
was calculated by the Schmidt graphical method Do] for run 17, selected 
as representative. The plot was based upon a thermal diffusivity of 
0.004 ft2/hr and an average surface conductance of 14 Btu/(hrft2°R) from 
the wninsulated tank results. The results of the graphical solution are 


shown in Fig.e l4a. The capacitance at the end of each time interval was 


calculated from the relationship, 
Q. = e AT (24) 
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Where Q, is the heat stored in redwood as calculated from the internal 
temperature differences and AT is the difference between the tank air 
and ambient temperatures. The resulting increase in insulation C,* with 
time is shown as curve (1) on Fig. 14b. Next, values of C,* for the 
liner plus the capacitors were caluulated from equation (8a) at each time 
interval using experimental values of T* and M*, These are shown as 
curve (2) on Fig. 14b, The increasing values of C,* for cae liner alone 
were then subtracted from the experimental values of C,* for the liner 
plus the capacitors, resulting in curve (3). This curve is essentially 
constant at C,* = 10.0 and thus is in agreement with the measured value 
for the added capacitors independent of the insulation. 

Similar treatment of the other charging runs would be expected to 
yield the same results, Scatter of the experimental data may be attri- 
buted to variances in conductance, the smaller flow rates showing less 
effects of heat transfer to the insulation, It is seen then that equation 
(8) gives good results for the case where the capacitance temperature 
follews that of the gas. For the more general application where external 
heat transfer (NTU, 0) enters, the determination of natural convection 
conductances for the external tank and ambient air should present no 
problem, 

7.3.3. Blowdown of Insulated Tank with Added Capacitors 

The mass of air in the tank prior to blowdown was M)= 9.0 lbs., 
resulting in C,* «© 1.44 from equation (5). This value is much lower than 
the infinite capacitance value of C,¥e 6.5. As with charging, the 
values of NIU listed in Table 10 have been estimated from the results 


listed in Table 6, p. 25, for the uninsulated tank. 
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TABLE 10 


Blowdown Conductances 


Added Capacitors in Insulated Tank 


fio. (Qbsyhr tetas ee E ite CS 
19 2551 2.9 8a6 1.44 
20 356 3.6 Ga 1.44 
21 497 4d tac 1.44 


Since the NTU values are each above the minimum value of NTU >7.5 
recommended by Reynolds [2] and since the capacitance temperature closely 
followed that of the air, the solutions for negligible inside resistance 


. 
should apply. If NTU = 0 and T, =1, equation (9) reduces to: 


ae_ (mt ace ot 
a) fio (9a) 


The solution of equation (9a) for C, « 1.44 is shown graphically 
on Fige 11. The curve follows the experimental points quite well for 
the higher flow rates with a gradual departure toward the end of the 
rune As in charging, the Byeten cartoaeknes ratio, oo would be ex- 
pected to increase with time due to heat transfer with the insulation. 
Such a change would be characterized by lower values of T with time. 
In blowdown, however, the initial capacitance of the air (My cy) is 
much larger than for charging. Thus slight increases in C, with time 
would have considerably less effect for blowdown than for charging. 
Relatively fast blowdown runs were made in order to minimize the heat 
transfer with the insulation and, as shown on Fig. 11, the fastest run 
(run 21) follows best the, theoretical curve. 


It would be expected that the slower runs would lie above this 
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curve in a region of increased C,* (lower T*),. That this is not the 
case is shown by run 19. Since the heat transfer conductances are much 
smaller for blowdown than for charging, the insulation surface tempera- 
ture does not follow that of the air and as a result the insulation acts 
much as an isothermal sink and source. Thus the slower runs appear 
similar to run 12 shown on Fig. 9 with the same behavior as discussed 
previously. All runs, however, exhibit the effect of heat transfer from 
the higher temperature insulation to the lower temperature air. 

Despite the anomalies introduced by the insulation, it would appear 
that the blowdown solution for negligible internal resistance yields 
reliable predictions for system behavior where the capacitance temperature 
follows that of the gas. As with charging, the introduction of external 
heat transfer should present no problem in the determination of heat 
transfer conductances. 

7.4. Mixing Considerations 

As stated previously, the derivations of the solutions evaluated 
in this report are based upon the assumption that the gas temperature is 
uniform throughout the tank at any given time. This assumes that during 
charging the incoming gas expands and diffuses instantaneously throughout 
the gas in the receiver. Therefore, it was considered expedient to 
perform an experiment during a representative charging run to determine 
the validity of this assumption. The test apparatus was so arranged 
(Fig. 2) that the four air temperature thermocouples could be monitored 
individually during a run. Run 4, selected as representative, was 
repeated under identical conditions, each time recording the output of 
a different thermocouple. The temperatures so obtained revealed no 
variance greater than 0.5 OF. This was indeed the expected result 


considering the mixing potential of the jet flow described earlier. 
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This is, however, contrary to findings briefly referred to by 
Reynolds [2] to the effect that considerable temperature gradients 
existed throughout a charging run. Although the geometry of the receiver 
was not described, such results might be attained with a tank where the 
charging jet cone diameter exceeds that of the tank. Such an arrange- 
ment might conceivably result in a piston effect causing temperature 


gradients between upper and lower layers of gas in the tank. 
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8. Gonclusions 

The conclusions to be drawn from this investigation are summarized 
as follows: 
(a) Excellent quantdtative agreement between the experimental results 
and the simplified analytical solutions of Reynolds [2] have been 
obtained for both charging and blowdown. 

_-(b) Constant mass charging directly through an orifice into the bottom 
of a receiver of sufficient diameter for the formation of a freely exe 
panding fluid jet results in forced convection heat transfer conduct- 
ances between the air and the tank walls. These conductances increase 
as the air temperature and pressure change with time and vary in magnitude 
over the interior surface. 

(c) During relatively slow constant mass flow blowdowns, the heat 
transfer conductances tend to be constant and a time average free con- 
vection conductance gives a satisfactory approximation for a large cap- 
acitance system. This method becomes less reliable as the flow rate 
increases and a more detailed analysis, related to the geometry of the 
system, would be required, 

“(d) Analytical expressions have been found to predict forced convection 
conductances during charging that closely approximate the experimental 
values. At relatively high mass flow rates, these expressions can be 
used with available stepwise solutions as a basis for preliminary design 
of high capacitance systems, At relatively low mass flow rates, the 
conductances are of sufficient magnitude to give essentially isothermal 
behavior. 

(e) Despite the high mass flow rates possible with a relatively large 


size air receiver, it is unlikely that an adiabatic charging process 
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would ever be attained within a reasonable charging time in the presence 
of a jet effect forced convection heat transfer mechanism. 

(f£) Closer approach to adiabatic behavior is attained during blowdown 
than during charging for the same constant flow rate. An approach limit 
is reached for blowdown, owen: ING even though the flow rate is 
increasing the heat transfer conductance is also increasing until the 
effect of the latter exceeds that of the former. 

(g) If the reservoir is of sufficient diameter for the formation of 

a fully expanding fluid jet during charging, there will be adequate 


mixing with no temperature gradient in the receiver. 
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9, Recommendations 

Tt is recommended that the general problem of charging and blowdown 
be analysed for solution by analog computer methods. Despite the 
problems encountered with the simplified computer setup evaluated in 
Appendix III, it is considered that the general problem of charging 


and blowdown is inherently suitable to treatment by such methods. 
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APPENDIX I 


Derivation of Charging and Blowdown Solutions 


1. Equations of State 


In all of the analyses to follow, the working fluid is presumed to 
be a perfect gas, with constant specific heats, having the following 
equations of state: PV = (R/m) T 

u= CL T 

Ay = cy At 

fe Cp Je 
(R/m) = 0, 7 Cy 


2. General Differential Equation for Blowdown and Charging 


a She 
dt 





An energy balance for the gas in the receiver gives: 
hea +h dM 
a dt 


An energy balance on the capacitance gives: 


ie 
pee 


Rate equations for the two heat transfer resistances give: 


ae ME — Te me Te -7T_ 
ban oe ; ak bi =~ R; 


The change in energy storage in the gas due to mass flow out 


and changes in specific energy is: 


_ d(Mu) _ JM 
4 at! =m ce+ug 
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Assuming the specific heat of the capacitance is constant, the change in 


energy storage in the capacitance is: 














dUc IT 
— Cc 
dx Cc. dt 
Combining the energy balances, rate equations, and the perfect gas equations 
of state: 
dT _fy_ dM y Te-T 
- CUNT 4 2 = 


Now since we. a » then; 


dT — dt dM ey eee ae dle ~ dle dM __, dle 
dt JM dz dM dX dM d& dM 


Thus inserting these relationships we may write: 


a / / ae 
M ha [(e-i) + OAM T T Ry Cw T. =O 








Helis |e! eee eee 
i dM Ge Fee: aa RQ Riou : 


These equations expressed in terms of dimensionless parameters become: 





ur M™ dT - [e-I)ur* + wry] T* + NTU (fis 





* 
and us” d Tc =e [wru.,+ NTULT, +g NTU T “45% MTUT <0 
GM” & G. G <0 «0 
’ 
These may be combined into a single equation ens as follows: 


«2 a ie’ ef a« 
ur Moe + ur {™ dur” 4 ur"(2-b)- -wru ~¥, A, (NTUFNTUs yf of 


dM" qT 





yf 
E e f (4-1) ur * (NTU + NTU, )#NTU wTU,- ur “(h-1) aT “- Mus C= 0 


x Bey 


Ad 





For a ie mass flow, w* = 1, Brae 





=" M" JT” 
da [eke wre BE (ura nrg) | 
* NT U 
at (nanan ei oe y= O 


Similarly, for charging at constant mass flow: 


2 


M gM? 
vA TU 
aa [NTU —- C,* san 








NTU NTU * 2 
— IY NIV = 
ae T., 


oO 


3. Adiabatic Charging and Blowdown 





An energy balance on the gas in the receiver gives: 


charging: a dM d(Mu)=Mdu +udM 
blowdown: f di d(Mu) = Mdu+u dM 


Combining the energy balance and the equation of State: 
charging: aap ae) dM MdT 
blowdown: (k-1) TdM MdT 


Expressing in terms of nondimensional variables and separating: 


{\ 


| 


{\ 


I\ 


k . M* 
charging: enor, = dM" 
Cie. M* 


blowdown: [C4 a - / (h- ') du 


45 


Integrating from the initial values to conditions at a later time: 


are 
Charging: “=e wi _ ft (6) 
Blowdown: I = hee (4-1) (7) 


4. Charging and Blowdown with Negligible Inside Resistance at Constamt 


Mass Flow 





An energy balance on the gas and capacitance combined: 


rging: dM _ = JT dU. dM 
Charging ae 4 C. SC tM Se + u 


dt 
lowdown: dM = au du JM 
Meee” Wa ~ Cae aie Oe 


> ane 
Charging: hao ts Blowdown: 5 ao~ a ae 


Since, for charging: & j= IM/Yt and for Blowdown: wo = _~ dM [gt 


Then, for charging: gt gg gl Teno oT 


i) 


meee, GT. IT SMe ae 
for blowdown AS = =m FE > Wo gM 


Combining this relationship with the energy equation, rate equation, 


and the perfect gas equation of state, yields the following differential 


equation: 
ne: M+ Se) th Soc re ee EO Ea 
Cha i e | te —ee a pt ES — 
a ( Cy M Rw CY %S 





T 
O 


d 
low : Ce aa Tort __ 
eS (M+ Ze) dM + ee a (k-1)T 
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Expressing in terms of the dimensionless variables and were 
as iis 
& x 
TAG*NTUL)-#T,"-NTUG TS 1 M*+c,8 
x 
wedn: a - a 
(A-1+ NTU,)T - To NTU wo é riers 


Integrating from the initial conditions to conditions at a later time: 


‘ e+e eg 
Charging: [| ee (ltNTUo - &T," - NTU Mow) Cot t+M* tT 4NTUsTe 4 
1+ NTU 2 


Blowdown: . Moris h-l4+nTTWg 
Fi eal ed AY Pn Ba, 0) + To NTV@ (9) 


(2 -| + NTU, ) 
5. Charging and Blowdown at Constant Mass Flow with Heat Transfer 


To or From an Isothermal Source o 





An energy balance on the gas gives: 


Charging: Ah, 4 7 > dtu) 
du dM 
ere at 


A rate equation for the heat transfer resistance Ry gives: 


— 


Blowdown: Ei - - = 


Charging: Kei = T-Te Blowdown: fea —_ ies 
Ke = 


Since: for charging, Ww = dM/ dé and for blowdown, Weir dM/ at 


= for charging: cel == dT. dM =1, dT 
dt aM at aM 
for blowdown: eee dT : lM a ran 
d- adm at ° “dM 
Combination of this relationship with the energy balance, rate equation, 


and perfect gas equation of state, gives: 


4 


a7 


Changing: dT —_ 
anging aa eo Pa te 


dT | 
Blowdown: = dei ae 
° M nT |Z I+ s 


Pee se _ Te _ 
cy W% |v uy 


Ki Cy wm 
Expressing in terms of nondimensional variables and separating: 
aes 


*& 
Charging: J: dT 


mM * 
' 


, y a dM) 
T* (4+ NTU)—- RT*~— NTU Te" 
T* 
Blowdown: Fi JT“ 
é T.* NTU -Ck-! 


Te 
A Cy Us, 





M* 
MF 
eee A 
+NTU) T* OM a 
! M4 
Integrating from the initial values to conditions at a later time: 
x 
Charging: T= 


-(\4 Na 
AT“ + NTU TE -[RT "= I-NTU +NTU TIME ag) 
L+NTU 
R+NTU * 
Berea = C h-! + NTU-Te NTU M* + Te NTU | 


k—| + NTU 
Or, integrating over the time interval At: 


: (11) 
T*+AaT* F riawe 
, dT i J 4 
Charging: — = Oe 
Te T* Cl4NTU)— 2#T,-NTU Te* Ar Mv * 
Blowdown : ew ag Ss: ‘ Mt+4 Mt. 
dT + | dM" 
& oe =a 
x TENTU -(k-1+NTU) T* ary O 
Yields the stepwise equations: 
D Sais NTU Tt k M*+4M* ve 
ae aa iboats = Jo | ii ae 
Charging: AT =| 1 +NTU q it ( Me (12) 
Blowdown: % a 
t_ fre NTU | | (MAM \* oe) 
an “fr k-\+ NTU 7] (Ge) a 
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Photograph of Test Apparatus 
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a. Top View Showing \Wiring Details 
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APPENDIX IIT 


Solution of the General Differential Equation for 


Constant Mass Flow Blowdown by Analog Computer Methods 
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APPENDIX III 


Solution of the General Differential Equation for Constant Mass Flow 
Blowdown by Analog Computer Methods 
1, Derivation of the Machine Equation 
Equation (2a) (Appendix I) may be rearranged in the case of 


NTU 25 = Q as follows: 


o ~~ * t 
dJ°T*, & dT a ie er O 





ae oe 2b 
dM* ~ M* gmt ° dM*” Oe | 
The constants xo, ® are defined in terms of the system parameters 
as follows: 
AH = 2-R-NTU (25) 
ay = NEE (26) 
Co 
= NT 
a= (4-1) ale (27) 
ve) 
Equation (2b) is readily adaptable for solution by analog computer. 
Adopting the standard nomenclature of Korn Gal: the variables of 
equation (2b) may be expressed in terms of computer voltages and 
functions as follows: 
— « 
7 =@-7 (28) 
eee a 
M =a M (29) 
M 
d 2a daa p (30) 
dM* © dt t 


In the above transfer functions, a, and ay are scale factors relating 
the machine voltages T and M to the variables T* and M* and a. is the 
time-scale factor relating differentiation with respect to the variable 
M* to the standard operator p denoting differentiation with respect to 


machine time t. 
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Substituting the transfer functions in equation (2b) results in the 


following machine equation for the general case of constant mass flow 
b lowdown : 


~~ —— = 2G! _ oa Sees 
Cae 2 7 “2 (a PT +27 \4o pT (2c) 


The initial conditions to be employed in the integration of equation (2b) 
are given by Reynolds [2] as follows: 


* * A x 
ALM, = 1: To = 1 and eae st 
dM /, 
2. 


(31) 
Block Diagram for Analeg Computer Setup 


Machine equation (2c) combined with the initial conditions is 


represented in block diagram form for solution by d-c analog computer 
as follows: 


= pile 


a.p T [= 


(ApT + 22 )+tpT 


ry 


= & 


ie 





Borgia: 
CONST. 


Figure 15 


Analog Computer Setup for Solution of General Equation 
for Constant Flow Blowdown 
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5. Evaluation of Analog Computer Solutions 

The analog computer block diagram of Fig. 15, suitably time and 
magnitude scaled, was set up on a Boeing Model 7079 Analog Computer with 
the output to a Reeve servo plotting board. Itwas originally planned 
to compare the computer blowdown solutions with those obtained from the 
simplified blowdown equations (Chapter 5) best fitting the experimental 
data. The solutions, although capable of qualitatively reproducing the 
general solutions shown in Fig.! of reference [2], were not sensitive 
enough to changes in the system parameters to prove reliable. 

Solutions for both the case of infinite capacitance and the case of 
negligible internal resistance were attempted with the computer. In the 
infinite capacitance case in which @ and & both approach zero (equations 
26 and 27), the computer plot showed a barely discernable change over a 


range of several integer values of NTU in equation (25). The same held 


true for the case of negligible internal resistance. Here, since C. has 

a finite value, changes in NTU affect BS and 6 as well as &. Once again, 
however, quite large changes in the value of NTU had little effect on the 
displacement of the computer solution. 

This lack of sensitivity is in part due to the fact that the exper- 
imental data lies, in all cases, in a narrow range of T* and M close to 
the origin of the graphical plot. In this range, Fig. 1 of reference [2] 
shows that there is little change in the trend of the general solution 
despite large differences in the magnitude of the system parameters. 
Attempts to enlarge the scale resulted in an unreliable performance of the 
plotting equipment. This program, however, could be of value where lower 
values of T and mM were of interest, since in this region the general 


solution diverges sufficiently to reflect smaller changes in NTU. 
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APPENDIX IV 


Experimental Data 
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Table 11 


Charging Run Data- Uninsulated Tank 


= 





Orifice: 1/3 in. dia, 


run 1 
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Co” = 300 Ty = 0.990 

1 sh P ; 
Sec) (OR ) (psia) on 
0 Be | ie 1.000 
5 5)13 1 pee ieee 
wie pial 19.5 Ls027 
20 58} 265 12029 
30 552 2509 1.029 
ho 552 BoA 1.029 
50 552 Boab 1.029 
60 552 1269 1.029 
30 552 51.28 1.029 
100 552 Coen 1.029 
20 552 696 1.029 
140 552 ioe 1.029 
160 552 665 1.029 
130 poe 9.5 £029 
200 52 101.8 1.029 

Run 2 
Orifice: 3/16 in. dia. 

C5 = 300 TI = 9.990 

t zi ak 
(sec) (2) (osia) oe 
0 Dea sys: 1.000 
10 562 26.7 rely, 
20 561 Sieg Le Ol); 
30 560 19.6 1.012 
ho 559 6067 1.0);2 
50 59 ELSO 1.012 
60 959 3006 1.01 
70 559 3320 lpioveal 
30 559 9567 1.01 
90 559 10220 neal 
100 559 106.3 Peon? 


W fe) 


00 160 
Neale ihaaUh 
1.32 e29 
1.65 eee 
1.95 1.90 
Ba2% O22 
2259 2260 
ee 2.32 
G50 300 
210 3093 
1.70 11.66 
5 30 Cae 
5.35 5260 
6.33 6220 
Ceo 6.68 
= 330 lb/hr 

Le = Osg22 

SF 
POO 1200 
50 142 
2.56 Zac 
Beo> 2uan 
1.10 3.93 
ersO leo: 
BS 5 etd 
6.00 ae 
Ca / Gece 
6670 663 
6.90 Oeoe 
1.00 Geog 
(Spe 6e02 





Orifice: L/h ivgeduas 


Cg = 300 

+t : 
(sec) (2 

0 533 

5 566 
10 575 
15 572 
20 569 
25 563 
30 563 
aD 567 
ho 57 
NS 567 
50 567 
55 pier) 


Run 3 
RS Pe 


TS = 0.990 
Vee : 
(psia) Tt 
Viet L000 
206i ieOoO 
5 lie 15 63 
2. 1.061 
S25 05° 
OAE 1.055 
(eae Nyon 
By aE Leos 
oe 1.6. 3 
13 tt 1.053 
i Bre 1.953 
Run ly 


Table in. 
(continued) 


@Pificé: 5/16 in, dia, 


C5 = 300 

iE ei 
(sec } (it) 
0 Dee 
265 56h 
560 pyin 
ie 576 
10.0 Bey 
12.5 Bier 
15.0 ie 
17.5 572 
200 Sy al 
25.0 569 
80.0 563 
Be.0 563 
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Wo = L1U:0 1b/hr 


ait 
1 
ae 


1.2900 
1.060 
1.079 
i638 
1.00) 
LebL 
IPS OE 
107 
Lire 
8 68 
1.069 
1.069 


Te = 02995 


1.00 
oO 
2.00 
22 
2693 
31:0 
Bel 
ho 
4.90 
pre) 
66/9 
teo5 


+ 


1.00 
1.32 
2.15 
a. 72 
Biss 
Soe) S 
eho 
Pros 
5.60 
6510 
6.6h 
tele 


1.00 
Lela 
Voss 
bee 
2ole 


3600 
ode 
ey 
5 oe 
O.38 
(eaus 





Table 12 


Blowdowm tun Data - Uninsulated ‘’ank 





Run 5 
wifice: 3/16 in, dia. Wo = 105 lb/hr 
o* = h1.o : Te = 06955 
t P 9 . 
(sec) (2) 2) (psia) a P at 
0 560 pe) 1.000 1.000 1.900 
1o 550 109.5 5Ou? 0955 0973 
20 5h3 195.5 0970 0920 0943 
ho 56 9506 Rood. ° 360 6906 
60 52u 9265 0936 0306 e362 
80 517 Otes 0921; Pryce 23 
100 SU; 81.9 0916 afl Fito 
120 503 ee 0903 670 o (8 
Lo 503 (Bes 0099 030 e (02 
160 502 65a 0396 59h 2663 
130 500 Gee e093 0550 Me llité 
200 193 53.8 2099 0513 0513 
220 96 ple Peis 70 0530 
2h0 Hoh hoe oooh 120 sit> 
ratte ace 
Orifice: 1/l, in, dia. Wy = 166 lb/hr 
Co = hil.o Te = 00955 
t 7 P i "i ” 
(sec) eo (psia) ae aa a 
O 560 11360 16000 1.000 1.000 
5 pee toe 2 09.303 0970 0932 
10 jys 107.9 0930 © 90 050 
al 543 Lh .2 ©9710 2916 095 
20 539 19,0 096); 0930 092 
30 Se: 96.1 0950 ye) es 2390 
ho 526 91.0 0936 gS00 w055 
59 515 3665 0926 e 760 e 321 
60 513 11,8 els e120 «(a0 
30 503 7300 0593 26)12 o fly 
100 196 63.3 bse 0566 2639 
120 1,90 Soe Se ol 95 0566 
140 135 Te? 0005 '20 ~500 





Table 12 


(continued) 
mun { 
OQrifice: 5/16 sigma. Wo = 29 (ec e/ir 
ot aE P. = On900 
t I P on ae, 
(sec) Age) (psia a pe 
0 56) 11),.8 1.000 12000 
5 55k 1090 0933 0950 
LO 5u7 Loe 970 °910 
15 Sho 100.2 0953 sth 
20 534 96.1 Pe) its) . 30 
25 530 920 090 AIO 
30 525 3 3olt Ret ae. 
BO) 520 28 0923 ° 70 
ho 516 due eSl5 ° 108 
WS 512 175 0907 0676 
50 507 (SP 0906 2 6 
55 Soh (eAe 0393 2516 
60 500 672 Fool 05J6 
65 97 6.0 e681 0550 
70 ho 60.6 03716 Fo tle: 
Run 3 
Orifice: 3/3 in. dia. Wo = 375 lb/hr 
Co = 11.0 TM 0.975 
ie fs P . a 
(sec ) i) (psia) 3c8 a 
O 53 Get 1.000 1.000 
5 512 110.3 09 39 0950 
10 53h Lohely Beas 900 
15 bak Cat AoE 0 359 
20 sl 0362 0950 203 
25 513 O62 0935 e (60 
30 507 3303 092k; 0717 
ho 97 Thed 2906 2 Oli 
69 77 sane “08 ls30 
65 475 51.3 ENG ol y2 
70 7h Heme 2 36h, oh Oly 


ve 





~“-~ 

1 
BER BC AMeUnNneKo 3 ct 

Sd 


Table 13 


Charging Run Data- Insulated Tank 


Run 


\O 





Qrifice: 5/16 in. dia. 


Ty = 9.998 
ih P ‘ 
(92) (psia) ine 
532 ao 2000 
565 2360 aleyele! 
532 32eu 1.09); 
592 hOels Leds 
595 lr ied 12126 
60h Eye jeleas 
697 66.0 ge ral 
610 (as Blais 
612 J3eL 1.9 
613 91.3 eee 
61h LOLs nls 
Run 10 


Qrifice: 3/3 in, dia, 


a P 

(M2) (psia) ey 

523 Lies 1.000 
545 1934 14030 
562 2301 2060 
573 OVer ues 
533 3h.» 3 Leis 
59h Todt ee 
601 Me 5 1.136 
605 Con ey 
693 (303 Path 
611 35.3 La 
612 96.6 10159 


i 


oi 


I 


e e ¢ @ @ 
Cong Gs Ole sors SN 


On ONUIAN Emo lw NN FE Ee 
® 
Cah SO ON a) 0) (CO FC 


1200 
122 
1256 
oS 
PY 
Pet5 
3620 
3.63 
hi. LO 
e95 
5.30 
OAs 


Wo = 1130 lb/hr 


Wo = 1320 lb/hr 





Table 13 
(continued ) 


hun 21 





rifice: 7/16 in, dia. Wo = 1960 lb/hr 


74 


T, = 1.00 

ft ab Pp Bs ae 
(sec) (Sf) (psia) _ vie 

9 530 aes 1.000 ney 8,6) 

zi 53k iipee 15032 Lees 

2 539 ES 1.065 Sh 

3 pou 3060 1.092 2.3 

yy 576 35.5 allie DEG 

5. Ca 12.6 eab2G Baal 

6 595 19 oh 1.133 3230 

7 602 eras Leet helo 

3. 693 61.6 1.156 1. 30 

9 612 (erga ios See: 
10 615 (HES 16166 Sane) 
11 617 ose 1,165 G.03 

Table Di 
lowdown .un Data=- Insulated Tank 
Nya bho ape 
Orifice: 1/) in. dia. Wo = 395 lb/hr 

t ci P ris iy 

0 599 ae 1.000 1.000 
20 555 ed ACES Bt 30 
ho 527 660 2 380 el 
50 B13 5Oeu o dk hie 
69 512 51k 53 ellis 
70 597 ood r) ah 3°20 
3) 502 39 03 0 339 BV) 
90 590 35% Be bie, 0307 
Le, 97 erat 2 230 e2nh 
120 96 Peo es, 0232 


1.900 
Fi sa) 
0/96 
0 {19 
0 hd 
so 5 
ace 
oliG2 
el 12 
0 308 
0337 
0306 
252 





Table 1) 


(contimed) 
mun 13 
wifices 3/3 in, dla. Wo = 31 lb/hr 

ie oh P 7 . e 
(sec) ie) sia ar ae = 

0 oe 1 1.000 1.000 1.000 

2a5 Sie lola 0935 e910 092); 

5 567 Soe 97 Bie af6l 
10 ES: (Ero 0930 (Oo 0705 
alles Sh noeo 0925 aoe) Pupils 
20 523 6002 090), od 0598 
30 508 eed 370 2lt09 et70 
ho hO6 Bere 3448 316 0375 

Run ly 
Orifice: 7/16 in, dia. i 360 lb/hr 

t vy P ie 2 % 
(sec (PR) (psia) pe a nS 

0 615 9363 1.000 1.000 1.090 

Beals 602 a7 20 AO ie e 338 091 

5 537 7326 095k 0796 0338 

feo 57h (Osc PS ee: ays (66 
12.5 554 OTe 0990 ABTS 02 
15 Shh SOag 2305 eo 0530 
eS 536 oe7 AO 0lt63 e508 
20 523 11.0 0859 ol15 033 
AAS 521 Bie! e346 ote olily3 
25 54 3L.0 s035 o 3h lite 
2705 508 culeal S006 035 2 301 
30 503 205C 2810 5293 oa58 
6265 590 26 Mohs 0273 5356 
35 97 2.8 2309 sal otal 


15 





sable 1) 
(continued ) 


Run 15 


Orifice: 1/25. Wo ~ 1100 lb/hr 


76 


T P J. : 5 

(sec) (Ox) (psia) a Pe 
0 616 Liss 1.000 12026 1.000 
2 69h Lipsy 0930 =i 0925 
y 592 Ohe7 2960 0 325 0358 
6 S31 35.0 sole ue 0706 
e 569 (ome 0922 NS 6 of it 
10 556 6iae 0902 0590 065) 
2 Oh7 59a/ Povey e520 0536 
UU; 537 Cie if 0372 C 0027 
16 530 Typ 0858 ‘ nd 
18 522 2.5 o8h16 0370 eet 
20 Sus: 38.3 0835 “Sie 2/400 
22 510 34.8 0828 303 0367 
2h, 505 3163 320 ecl3 e205 
26 h99 26a 2810 20 © 309 

Table 15 
Charging itun Data- Insulated Tank 
Wi Aaae apaclLtvors 
Run 16 
Orifice: 3/16 in, dia. Wo = LO7 lb/hr 
GC, = 10.0 Ty = 1.00 

% a P 2 2 

(sec) (PR) (psia) le Br. NM 
@) 530 1.3 1.090 1.00 1.90 
5 533 C220 1.016 50 13 
10 5h5 29.6 eae 2.00 1.9): 
15 552 Bi A6 TeOhO 2050 Poe 
29 556 Wiel 1.050 3200 285 
25 561 51.8 1.059 6250 Brot 
30 56h 5902 1.065 1.00 3279 
B5 567 67 el 1.070 We D 629 
ho 570 The6 Alaa SAC eal 
LS 571 8252 1.075 Beo5 5 el7 
50 572 3965 1.089 6.05 5260 
60 575 10320 1,085 6.95 6eh10 





Table 15 
(continued ) 


Run 17 


Orifice: L7i tan scae 


CS ~ Oz Ty _ 1.00 
7 P : , 
(sec)  (°) (psig) oat 
0 527 1.8 1.000 1.00 1.00 
5 5h. 2607 1.026 1.80 1S 
10 550 3365 V0 260 2019 
15 S56 Saal 1.056 Sous Bre 
20 562 6346 1.065 1.30 Oly 
25 566 7720 1.073 5620 Shy 
30 570 91.0 1.081 6615 5.68 
Bb ie 10326 1.088 7000 | 60h3 
ho S77 L369 1.095 HAL 70h 
Run 18 





Opiftied: 5/16 in, dia. w, = 110) lb/hr 


Cy = 10.0 T” = 1,00 
1 

t 3 P : . 
(sec) (PR) (psia) T oe Mu 
0 531 1.8 1,000 1500 1.00 
205 538 2he0 1.017 EY. 1.60 
S Sh7 Boao 1.032 2.25 Dale 
eo Bes 1269 1.05 2.90 2.77 
560 52 elt IOs San 2oo5 
12.5 565 622 1.06) 4220 3695 
15 563 (OA 1.071 4.35 053 
es 571 O16 1.076 Eee Sys 
20 57h 91.6 1.081 6.20 5.73 
22.5 576 192,2 1.086 6.90 6,36 
30 578 112.0 1.090 Te00 6.94 


al 





Table 16 


Blowdown Data - Insulated J'ank with 
Added Capacitors 


Run 19 


Urifice: 5/16 in, dlae Wo - 255 lb/hr 


78 


Cy = Alen 
ir Is ot Ate 
Ce) esta) EOP 
Sis 111.5 1.000 1.000 
567 10367 0990 0930 
561 98.0 2980 2880 
556 91.) 0971 2820 
552 O7e0 964 e180 
5,8 8023 0955 areD 
ohh (GAs °950 0105 
539 (ee o9h0 Pioys 
534 6S55 0931 0570 
eoel 6022 0926 2540 
529 Boe! 0922 °500 
526 51.3 0918 60 
Sah 1543 09 Ly Any 
Run 20 
Qrifice: 3/8 in, dia. Wo ~- 356 1b/hr 
Co - Leb 
T P 

() (psa) Ba 
577 112.6 1.000 1.000 
570 LOhe5 0989 © 930 
565 646 0919 2360 
553 88.7 0968 2790 
553 820 0959 #130 
53 102 9 ACHE: 
538 Coa 0934 0580 
535 5926 Piied| 0530 
Soi 5260 0920 0470 
527 6.6 Feabn B15 





Table 16 


19 


(continued) 
Run 21 
Orifice: 7/16 in. die. Wo - 497 lb/hr 
oF ao Lely: 
+ P e = 
(sec) (2) (pata) TES 
on eis Taal 1.000 1.000 
205 568 105 elt 0993 095 
5 565 100.3 0986 0900 
as 560 Ded 0973 Ee 
10 Sil 90e1 0972 2310 
TAS 55h 8505 0966 » 168 
25 559 Blige 0960 ° 730 
765 55 7626 B52 2638 
20 5he2 720k 096 2550 
22.5 538 6729 090 2610 
2 536 63.8 0935 ARG 
Bee Ds 596 0930 0536 
30 530 oe 0925 0500 
32.5 527 51.2 0920 1:60 
Bo 525 7.2 0916 ellen 























‘ 
em 


